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Abstract. We consider semisimple triangular operators acting in the sym- 
metric component of the group algebra over the weight lattice of a root sys- 
tem. We present a determinantal formula for the eigenbasis of such triangular 
operators. This determinantal formula gives rise to an explicit construction 
of the Macdonald polynomials and of the Heckman-Opdam generalized Jacobi 
polynomials. 
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1. Introduction 

The main objective of this work concerns the explicit computation of families 
of orthogonal polynomials associated with root systems. Key examples of the 
families under consideration are the Macdonald polynomials |M1I IM2I IM3| and 
the Heckman-Opdam generalized Jacobi polynomials HS , 0|. The origin of the 
Heckman-Opdam polynomials lies in the harmonic analysis of simple Lie groups, 
where they appear (for special parameter values) as zonal spherical functions on 
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compact symmetric spaces |HSI IHej . Other important applications of these poly- 
nomials arise in mathematical physics, where they are used to express the eigen- 
functions of the quantum Calogero-Sutherland one-dimensional many-body sys- 
tems |S1I IS21 IOP1 IHS| . The Macdonald polynomials have similar applications: 
they appear as zonal spherical functions on compact quantum symmetric spaces 
[NllSu) . and they are used to express the eigenfunctions of Ruijsenaars' (<7-)difference 
Calogero-Sutherland systems |Rl ID1) . Depending on the specific application of in- 
terest, our work may thus be viewed as providing an explicit construction for the 
zonal spherical functions on compact (quantum) symmetric spaces or for the eigen- 
functions of the (difference) Calogero-Sutherland type quantum many-body models. 

The usual definition of the Heckman-Opdam and Macdonald polynomials in- 
volves a Gram-Schmidt type orthogonalization of the monomial basis with respect 
to a generalized Haar measure HS ; [01 IM2I |M3 . This definition, although most 
appropriate from a theoretical point of view, is not very adequate for the explicit 
computation of the polynomials in question. The main result of this paper is a 
determinantal formula for the Heckman-Opdam and Macdonald polynomials that 
gives rise to an efficient recursive procedure from which their expansion in the 
monomial basis can be constructed explicitly. For the type A root systems the 
Heckman-Opdam polynomials reduce (in essence) to Jack's polynomials (HH IM1| 
and the Macdonald polynomials reduce to Macdonald's symmetric functions |Mlj . 
In this case the determinantal construction of the polynomials under consideration 
was laid out in previous work by Lapointe, Lascoux, and Morse jLLMll [ELM2 . 
More specifically, the results of the present paper constitute a generalization of the 
methods of Refs. [LLM1I ILLM2j to the case of arbitrary root systems. For the 
Heckman-Opdam families we consider general (not necessarily reduced) root sys- 
tems and general values of the root multiplicity parameters. For the Macdonald 
families, however, we restrict for technical reasons to those (reduced) root systems 
for which the dual root system R v has a minuscule weight (thus including the types 
An, Bn, Cat, Dn, Eq and £7 while excluding the types BC n , Eg, F4 and G^)- 

The paper is organized as follows. Section [21 introduces the concept of a trian- 
gular operator in the Weyl-group invariant component of the group algebra over 
the weight lattice of a root system. In Section [3| we present a method for di- 
agonalizing such triangular operators by means of a determinantal formula. The 
Heckman-Opdam and Macdonald polynomials are defined in Section^] We employ 
the determinantal formula from Section[31to build explicit expressions for the mono- 
mial expansions of these polynomials in Section^ (Heckman-Opdam) and Sections 
El (Macdonald) , respectively. For completeness, some technicalities concerning 
the explicit evaluation of the determinant of a Hessenberg matrix are recalled in 
Appendix 1X1 at the end of the paper. To facilitate explicit computations, we have 
furthermore included a useful formula for the calculation of the orders of stabilizer 
subgroups of the Weyl group in Appendix iBl 

2. Triangular Operators in the Symmetrized Group Algebra 

In this section we define the concept of a triangular operator in the Weyl-group 
invariant component of the group algebra over the weight lattice of a root system. 
For preliminaries on root systems the reader is referred to e.g. Refs. |Bl IHuj . 

Let E, (•, •) be a real Euclidean space spanned by an irreducible root system R 
with Weyl group W. We write Q and Q + for the root lattice and its nonnegative 



DETERMINANTAL CONSTRUCTION OF POLYNOMIALS 



3 



semigroup generated by the positive roots R + 

Q = Span z (i?), Q+ = Span N (i2+). (2.1) 
The weight lattice V and the cone of dominant weights V + are given by 

V = {X e E | (A, a v ) 6 Z, Va £ R} (2.2) 

and 

?+ = {Ae_E|(A,a v )eN, Vaefi+}, (2.3) 
where a v —2a/ {a, a). The weight lattice is endowed with the natural partial order 

A h ii <^=> A - n e Q+. (2.4) 

Let Q v denote the dual root lattice generated by the dual root system R y = 
{a v | a £ R}. The group algebra over the weight lattice R[P] is the algebra gen- 
erated by the formal exponentials e A , A e V subject to the multiplication relation 
e A e M = e A+M . This algebra can be realized explicitly as the algebra A of (Fourier) 
polynomials on the torus T = E/(2ttQ x/ ) through the identification 

(with i£T). Symmetrization with respect to the action of the Wcyl group produces 
the basis of monomial symmetric functions {m\}xev+ f° r the space A w of Weyl- 
group invariant polynomials on T, where 

m x = eM ' A e P ' ( 2 - 6 ) 

with W(X) denoting the orbit of A with respect to the action of the Weyl group. 

We write A^ for the finite-dimensional highest weight subspace of A w with 
highest weight A e V + , i.e., Af = Span{m^}^ eP+;([t ^ A . 

Definition. A linear operator D : A w — > A w is called triangular if D(Af) C 
for all AeP+. 

3. DETERMINANTAL DlAGONALIZATION 

The triangularity of a linear operator D in .4 147 reduces its eigenvalue problem 
to a finite-dimensional one. In this section we diagonalize the triangular operators 
through a determinantal representation of the cigenfunctions. 

Let D be a triangular operator and let {s\}\ €V + be a second basis of A w that 
is related to the monomial basis by a unitriangular transformation: 

™ A = X! a ^ s P> a AA = 1 (3.1) 

(A e The triangularity implies that the expansion of Dm\ in the basis 

{s\}\ e -p+ is of the form 

Dm x = ^ &a m s m , &aa = ex, (3.2) 

(USP+./^A 

with the diagonal matrix elements ex, A e "P + being precisely the eigenvalues of D. 
Definition. The triangular operator D is called regular if ^ e\ when /j, -< X. 
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For a regular triangular operator the eigenvalues e x , A G V + are semisimple. Let 
{p\}\ev+ be a corresponding basis of eigenfunctions diagonalizing D. Clearly, p x 
has a monomial expansion of the form 

P\= C V m A" C AA = 1, (3.3) 

where we have normalized such that p x is monic. The following theorem provides 
an explicit determinantal formula for p x , given the action of D on m A expressed 
in the basis s x , i.e., given the expansion coefficients a Xfl and b Xfl in Eqs. and 

Theorem 3.1 (Determinantal Formula). Let D be a regular triangular operator in 
A w whose action on the monomial symmetric functions is given by Eqs. and 
(|3.2(l . Then the monic basis {p\}\ev+ of A w diagonalizing D, in the sense that 

D Px = e xPx , VAg7>+, 

is given explicitly by the (lower) Hessenberg determinant 



e A< 2 ) - e A(") . . . 



: •• 

d A (n-i) A (2) • • • e A( „-i) - e A („) 

Here A' 1 ' < A' 2 ' < • • ■ < A'™ -1 ) < X^ — X denotes any linear ordering of the dom- 
inant weights {p G V + \ p ^ A} refining the natural order 1)2. 4|) . the normalization 
is determined by 

and the matrix elements d x u) X (k) (n > j > k > 1) read 

d x u)\m —b x u) X (k) — e\a X (j) X ( k ) . 

Proof. Expansion of the determinant with respect to the first column produces a 
linear combination of monomials in the highest weight space A^ . The coefficient of 
the leading monomial m x is given by (— l)™ -1 times the product of the elements on 
the super-diagonal, which are nonzero by the regularity condition on D. Division by 
£ x thus gives rise to a monic polynomial. It remains to show that this polynomial is 
an eigenfunction of D with eigenvalue e A . To this end one observes that the action 
of (D — e x ) on the determinant affects only its first column. Indeed, we get — upon 



m X (i) e A( i) - e A (n) 
rn A <2) d A (2) A <i) 



PA 



m X ( n -i) G? A ( n -i) A (i) 

m A (n) <f A (n) A (l) 
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invoking the expansions 1(3. 1|1 and 1(3.2(1 — that 

■ ■ • 2ft=l d X (f)\(k)S X (k) + (e x u) - e\)s x u) 

d x u) X (i) 



(D - e x ) Px = — 



d x u) X u-i) 

£a<j) — 6 A 




(where, for typographical reasons, we have taken the transpose of our matrix). The 
latter determinant has a first row of the form s A (ijrow2+s A (2)row3+- ■ -+s A („-i)row n , 
and thus vanishes identically. □ 

As a corollary of the determinantal formula for p x , one arrives at a linear recur- 
rence relation encoding an efficient algorithm for the computation of the coefficients 
c Xfi entering the monomial expansion 1)3.3(1 . 

Corollary 3.2 (Linear Recurrence Relation). The monomial expansion of p x is of 
the form 

n 

Px = ^2c xx( i) m xW , 
1=1 



with c XX ( n ) = c xx = 1 and 
(Ki<n). 



1 - 



£ A — e A«-i) 



k=i 



Proof. Immediate from Theorem 13.11 and the Hessenberg determinant evaluation 
given by the lemma in Appendix 1X1 at the end of the paper. □ 

Moreover, by solving the recurrence relation we arrive at the following explicit 
expression for the coefficients c Xtl of the monomial expansion l|3.3|) . 

Corollary 3.3 (Explicit Monomial Expansion). The coefficients of the monomial 
expansion p x — Ya=i c aa<*) m \( l ) are given explicitly by 

d x Uo) X Ui)d x ui) X U2) • • • d^j,,-!) X ( jr) 



c AA( f ) 



E 



-Or <jr-l <---<jl<jo=n 
r— 1 ,n—£ 



(£A - £a«i)) ' ' ' ( e A - £\Ur)) 



with the convention that empty sums are equal to 1 (so c XX ( n ) — c xx = 1). 

Proof. In view of Corollary 13.21 it suffices to check that the stated expression for 
c XX (i) represents the (unique) solution to the linear recurrence relation of Corollary 
13.21 subject to the initial condition c XX ( n ) = 1. Firstly, the convention that empty 
sums are equal to 1 guarantees that the initial condition is satisfied. Secondly, 
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by isolating the last factor in the numerator and denominator, it is seen that for 
1< i < n 



-AA(«-i> 



E 



d\(j ) \(ii)dx(ii) ' ' ' d x u r _ 1 ) X ( jr - j 



. ( e A - EaOi) ) ' ' ' ( £ A — e A(^) ) 

r=l,...,n-l+l 



E I d\(h) X 

d\Uo)\Ui)d X (i 1 ) X U2) ■ ■ ■ d x (j r - 



ex - e AC /-D fc= , 



2'A l 



y 

u ■ ^ ~? ^ ^ (ea - e A(3l) ) ■ • • (e A - e A o r _ l) ) 

r— 1=1 ,n — k 



J2 d 



ex - e AC «-i) fc=/ 



□ 



4. Orthogonal Polynomials 



In this section the Heckman-Opdam and Macdonald polynomials are defined. 

Let A(x) be a positive continuous function on the torus T/(2ttQ v ) that is in- 
variant with respect to the action of the Weyl group (i.e., A(w(x)) = A(x) for all 
w G W). We equip A w with an inner product structure associated to the weight 
function A 

(f,g}A = -r^-JjgAdx (f,geA w ), (4.1) 

where ~g denotes the complex conjugate of g and \W\ is the order of the Weyl group. 
Let {pa,a}ag , p+ be the basis of A w that is obtained from the monomial symmetric 
basis {mx}x£V+ through application of the Gram-Schmidt process with respect to 
the partial order y (12. 4|) . More specifically, by definition px,A is the polynomial of 
the form 

Px.A = c a^( A )"V> ca,a(A) = 1, (4.2a) 

fiev+.n^x 
with coefficients ca, m (A) such that 

(pa,a,"v)a =0 for fi eP + with fj, -< A. (4.2b) 

In the general, the basis {px.a}x£V+ is n °t orthogonal with respect to the inner 
product (•, -)a (|4.1|l . as the natural order y (|2.4|) is not linear (unless dim(E) = 
1). However, for two important special choices of the weight function A it has 
been shown that the above partial Gram-Schmidt process indeed does produce an 
orthogonal basis jHSl IHl lM2l lM8j : 

A HO = Y[{l-e a ) 9 ", (4.3a) 

aefl 

with g a > such that g w t a ) — g a , Vw € W, and 

(e a ;g), 

with (z; q)oo = YI™ =0 (1 — zq m ) and < q, t a < 1 such that t w ( a } = t a , Vw 6 W 
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Definition. The orthogonal polynomials p\,A, A £ V + associated to the weight 
functions A HO (|4.3al) and A M (|4.3b|) are called the Heckman-Opdam polynomials 
|HSI IU| and the Macdonald polynomials M2 , M3 , respectively. 

Remark (i). fn the above definition of the Heckman-Opdam polynomials we may 
allow for a root system R that is nonreduced. For the Macdonald polynomials, on 
the other hand, we always assume that R be reduced. 

Remark (ii). The weight functions A HO and A M are invariant with respect to the 
action of the Weyl group on the variable (i.e. A(w(x)) = A(x)) because of the W- 
invariance of the orbit parameters g a and t a . The VF-invariance of these parameters 
moreover implies that the values of g a and t a depend only on the length of the root 
a. 

5. Heckman-Opdam Polynomials 

In this section we apply the formalism of Section [3] to arrive at a determinantal 
construction of the Heckman-Opdam polynomials for arbitrary (not necessarily 
reduced) root systems. 

5.1. The hypergeometric differential operator. To a vector x £ E we asso- 
ciate the directional derivation d x in A, whose action on the exponential basis is 
given by 

d x e x = (X,x)e x (A e7>). (5.1) 

Definition ( jHSI 10] ). Let x±,. . . ,xn be an orthonormal basis of E. The second- 
order partial differential operator 

D = XX+ E (5 - 2) 

3=1 aeR+ v 7 
is called the hypergeometric differential operator associated to the root system R. 

Clearly the definition of D (|5.2|l does not depend on the particular choice for 
the orthonormal basis x\,... ,xn- It is known that the hypergeometric differential 
operator maps the space of invariants A w into itself and, furthermore, that the 
Heckman-Opdam polynomials form a basis of eigenfunctions on which the operator 
acts diagonally |HSI lO] . 

We will now compute the action of D H5.2J) on the basis of monomial symmetric 
functions. To this end some notation is needed. We denote by W\ C W the 
stabilizer subgroup of A £ V 

W\ = {w £ W | w(X) = A}, 

by p g £ E the weighted half-sum of the positive roots 

p s = \ E 9° a ' 

by r a : E —> E the orthogonal reflection in the hyperplane perpendicular toaei! 
through the origin 

r a {x) — x — (x, a J )a (x £ E), 

and by [p] the function that extracts the integral part of a nonnegative real number 
p through truncation. 
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Lemma 5.1 (Action of the Hypergeometric Differential Operator). The action of 
D (|5.2|l on m\, A G V + . is given by 

Dm X = ((X + Pg,\ + Pg) - (Pg,Pg))m X 

1 [<A,a v }/2] 

+ r^E(^< A > a > E \Wx- ea \\W a (X-£a)\m x ^ a ), 

w/iere W^ a C denotes the subgroup of order 2 generated by r a (so \W a (X — fo)| 
is equal to 1 i/^ 1 = (A, a v )/2 cm<i egita/ £o 2 otherwise). 

Proof. The computation of the expansion of Dm^ in the monomial basis hinges on 
the fundamental identity 

= (A,a)e*(l + e-»)( 1 _ — ) 

(A,a v }-1 

= (A,a)(e A + e A -< A ' QV > Q ) + 2 (A, a) ^ e A - £a 
= (X,a)(e X + e r «W) 

[<A,a v )/2] 

+ (A,a) \W a {X-ia)\(e x - ta +e r ^ x - ea) ). (5.3) 

Indeed, the following sequence of elementary manipulations reduces the computa- 
tion of the action of the first-order component of D l|5.2|l on ra\ to an application 
of identity (|S~3|) : 

E W_\-cA 9 « m > 

aeR+ V 7 



E mi— ^0 5 " E 

1 v-^ / v-^ / 1 + e- a 



E w ( E e —) d <* eX 



\W\\ 

1 1 -mew qg_r+ 

1 K^ v }/2] 

= TT7T7 E ^(g Q (A, a) e A + g Q (A, a) ]T |W°(A - Ml ^ 

' A ' wGW 1=1 
/ j [<A^ V >/2] 

E^(A,a)m A + WT E l M/ A-£a||^ tt (A-fe)|™ A _, c 
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Combined with the action of the second-order component of D Ij5.2(l on m\ 

N 

^2 d x 3 m A = (A, A) m\, 

this produces the formula of the lemma. □ 

It is a standard property of root systems that for any A G V + the integral convex 
hull V\ = {p G V | WQi) di A} is saturated, i.e., if p G V\ then p — £a G T x for 
every integer £ between and (/i, a v ) (extremal values included) |Hu) . Hence, it 
follows from Lemma l5.1l that the hypergeometric differential operator is triangular. 
To compute for p dominant the coefficient of m M in Dm\, it suffices to collect all 
terms in the lemma for which A-fa£ W(p). Notice in this connection that for a 
given a G R + the a-string A — a, A — 2a, . . . , A — [(A, a v )/2]a may hit the Weyl orbit 
of p at most once. Indeed, it is clear from expanding both sides of the equality 
||A-fa|| 2 = ||A-fo|| 2 that A — I' a G W(X - la)— with 1 < £',£ < [(A,a v }/2]- 
implics £' = I. 

We thus end up with the following explicit triangular matrix representation of 
the hypergeometric differential operator with respect to the monomial basis. 

Proposition 5.2 (Triangular Expansion). Let A G V + . We have that 

Dm\ = e\m\+ ^ bx^m^, 

with 

£ A = {X + Pg,X + Pg) - (Pg,Pg), 
i \Wu\ /-v \ f \ 

Am = YW\ ^ ,9a (A,a)nA M (a)• 
l 1 aelA.rf 

Here [A, p] denotes the subset of roots a G R + for which A-fa£ W(p) for some 
(unique) I G {1,2,..., [(A,a v )/2]}, and 



if |H| = ||P a (A)|| 

*/ IHI^II^(A)||' 



where P a = (Id + r a )/2 is the orthogonal projection onto the hyperplane perpendic- 
ular to a through the origin. (So for A — ia G W(p) we have that nx^{a) = 1 if 
£ = (A,a v }/2 and nxp,{a) = 2 otherwise.) 

The regularity of the hypergeometric differential operator is immediate from the 
following proposition. 

Proposition 5.3 (Monotonicity). For nonnegative parameters g a , the eigenvalues 
ex = (A + p g , A + p g ) — (pg, p g ) are strictly monotonous in A G V + , i.e., 

VA, p G V + : p -< X ==>■ < ex- 
Proof. Assume A, p dominant with p, -< A, and let = A — /x (so G Q + ). Then 

£a -e M = (v, v) + 2(p + p g ,f), 

which is positive in view of the fact that (v, v) > and (p,+ p g , v) > (since both 
and p ff lie in the closure of the dominant Weyl chamber, cf. remark below). □ 
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Remark. In the proof of Proposition 15.31 we used the fact that for nonnegative 
parameters g a the weighted half-sum p g lies in the closure of the dominant Weyl 
chamber {x G E \ (x,a) > 0, Va £ R + }. This follows from the fact that, for 
any root (3, the partial half-sum p((3) of positive roots with length \\(3\\ lies in the 
dominant cone V + : 

\\a\\=\\0\\ 11^-11=11/31 

where {aj}i<j<Ar denotes the basis of simple roots generating Q + and {^j}i<j<N 
is the corresponding dual basis of fundamental weights generating V + , such that 
a k) = ^j,k- Indeed, since any simple reflection r aj permutes the positive roots 



J • 



other than ctj |Huj . one has that 

r aj W)) 




if INI = ||/3||, 
otherwise. 



It thus follows, from working out both sides of the equality (r aj (p(j3)) 1 r aj {ctj)) = 
(p(/3),aj), that (p(/3),aj) is equal to 1 if ||ay|| = ||/?|| and is equal to otherwise. 

This entails that p(f3) = J2f=i (p(P)> a 'j) u} j = Y^jLi, ||a-||=||0ll G ^ + as claimed. 

5.2. Determinantal construction. It is known that the hypergeometric differ- 
ential operator D l|5.2[) is symmetric with respect to the inner product (•, -)a Q4.1[l . 
associated to the weight function A HO Q4.3a|) . Combined with the triangularity, 
this implies that the eigcnbasis {px}xev+ diagonalizing D is given by the Heckman- 
Opdam polynomials px a, A G T J+ , defined in Section0|through the (partial) Gram- 
Schmidt process |HSI 10] , Moreover, since D is regular (Proposition [OJ, and its 
triangular action on the monomial basis is known explicitly (Proposition 15 . 2]l . we 
can in fact construct this eigenbasis in closed form by means of the determinan- 
tal construction in Section [3] (with sx — m\, so ax^ = 1 if p, = A and ax^ = 
otherwise) . This gives rise to the following explicit representation of the Heckman- 
Opdam polynomials. 

Theorem 5.4 (Determinantal Construction). For A G V + , let 

P\ =mx+ ^ C V "V 

denote the (monic) Heckman-Opdam polynomial with parameters g a > 0. Then we 
have — upon setting for p, v G V + 

tp. = (/•* + Pg,P + Pg) ~ (Pg,Pg), 

with [p,,v] C R + and n flv (a) in accordance with the definition in Proposition \5.S\ — 
that: 

i) the polynomial px is represented explicitly by the determinantal formula in 
Theorem \S.l\ 

ii) the coefficients cx^ of its monomial expansion are generated by the linear 
recurrence in Corollary[ 
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iii) the expansion coefficients c\^ are given in closed form by the formula in 
Corollary \EM 

Given a concrete root system R, Theorem l5.4l turns into an efficient algorithm for 
the computation of the associated Heckman-Opdam polynomials. We will illustrate 
this below for the classical root systems. 

Remark (i). The orders of the stabilizers in Theorem l5.4l can be computed by means 
of the formula 

t-t (p, a v ) + 1 + i<5° 

<A : q v )=0 

where p — \ X)qG-R+ a ' anc ^ ^ a ~ •"■ ^ a ^ ^ anc ^ zer0 otherwise. This expression 
can be found e.g. in Ref. |M3I Section 12], where it appears as a special case of 
the norm formulas for the Macdonald polynomials. For the reader's convenience, 
we have included a short proof of this formula in Appendix FBI 

Remark (ii). It is clear from Proposition 15 . 21 that the matrix b\^, representing the 
action of the hypergeometric differential operator with respect to the monomial 
basis, is quite sparse. The same is therefore true for the matrix d M „ appearing 
in the determinantal formula for the Heckman-Opdam polynomials in Theorem 
15.41 This means in practice that the algorithm for generating the Heckman-Opdam 
polynomials with the aid of the determinantal construction turns out to be much 
faster than one would expect based just on the size of the matrices involved. 

Remark (iii) . It is immediate from Thcorcm l5 . 4l that the coefficients in the monomial 
expansion of the monic Heckman-Opdam polynomial p\ are of the form cx^ — 
P\ti(9a)/q\^{g a ), where px^ga) and <3V(5a) are polynomials in the parameters 
g a that have nonnegative integral coefficients (and with the denominators q\fj.{g a ) 
dividing the normalization factor £\ = Y\^ v + u^xi 6 ^ ~ £ m))- Recently, a much 
stronger positive-integrality result for these expansion coefficients was found by 
Sahi [5a] (see also jMU ESI for the case R = A N ). 

Remark (iv). It is clear from the proof of Proposition 15 .31 that the hypergeometric 
differential operator D 15.2f> is in fact regular for generic (complex) parameters g a 
such that (is, v) + 2(p + p g , v) =^ for all v G Q + and p 6 V + . Hence, the deter- 
minantal construction of the Heckman-Opdam polynomials in Theorem l5.4l as the 
eigenbasis for the hypergeometric differential operator, extends meromorphically to 
g a in the complex plane. 

Remark (v). It is known that the coefficients of the Heckman-Opdam polynomials 
can in principle be computed by means of (cumbersome) Frcudenthal type recur- 
rence relations HS . From this perspective, the determinantal construction of The- 
orem 15.41 thus provides the explicit solution to this Frcudenthal type recurrence. 
The recurrence in Part ii) of Theorem 15.41 - which arises as a particular case of 
the general recurrence scheme in Corollary 13.21 upon choosing for our triangular 
operator the hypergeometric differential operator — reads concretely 



n l ^ / ' r l 

,«) n X (k) X (e-i)(a) c AA(fc) . 



\W x w\ ^ 
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It may in fact be seen as a suitable symmetric reduction of the Freudenthal type 
recurrence relations, enabling their explicit solution in closed form via Corollarv l3.3l 
When g a — 1 , Va G R, our recurrence is closely related to the optimized Freudenthal 
recurrence scheme for the computation of weight multiplicities of characters of 
simple Lie groups due to Moody and Patera |MP| . 

5.3. Tables for the classical root systems. We will now provide tables of the 
matrix elements building the detcrminantal formulas for the Heckman-Opdam poly- 
nomials associated with the classical root systems. In each case, we will only list 
the minimum amount of information needed for constructing the matrix, viz., i. 
the cone of the dominant weights and its partial order, ii. the eigenvalues building 
the super-diagonal of the matrix, Hi. and the values of the lower-triangular matrix 
elements. For further data on the root systems of interest we refer to the tables in 
Bourbaki |B"]. 

It will be convenient to parameterize the dominant weights of the classical root 
systems in terms of iV-tuples 

A = (Ai, A 2 , . . . , X N ) (5.5a) 

of weakly decreasing (half-)integers (so Ai > A2 > • • • > Xn)- Often we think of 
these iV-tuples also as multi-sets of the form 

A = {Ai,A 2 ,...,Ajv}, (5.5b) 

where the parts Xj are listed from largest to smallest. For A = (Ai,...,Ajv) 
(= {Ai, . . . , Aat}), we define X £ = (Ai, . . . , Ajv-i, e|Ajv|) (= {Ai, . . . , X N -i, e\X N \}), 
with e 6 {1,-1}. We need the following two operations on our weakly decreasing 
iV-tuples: 

X\fi = {Ai,A2,...,Ajv}\{Mi>M2,---,Miv}, (5.6a) 
A M = (A+\Ai+ (/i+\A+) E ), (5.6b) 

where e = sign(Ajv) xsign(/^j\r). The first operation takes the difference of A and \i as 
multi-sets, i.e., taking into account the multiplicities of the parts. (By convention, 
we will list the parts of this difference again from large to small.) The second 
operation encodes — up to a possible sign — the symmetric difference of A + and \x + . 
For example: (5, 3, 2±, 1, 1) Q (4, 3, 3, 1, -1) = ({5, 2±}, {4, -3}). Finally, for future 



reference we furthermore introduce the operations 

|A| = Ai + A 2 H hAw, (5.7a) 

Vx(m) = \{j = 1, . . . ,N I Xj = m V X 3 = -m}\, (5.7b) 

A = (Ai, A 2 , . . . , Ayv-i, —Xn), (5.7c) 



producing, respectively, the sum of the parts, the number of parts with specified 
absolute value |m|, and the conjugate TV-tuple with the sign of the last part flipped. 

5.3.1. The case R = A^-i- For the type A root system the Heckman-Opdam 
polynomials reduce (in essence) to Jack polynomials [HH IMlj . The determinantal 
construction in Theorem 15.41 reproduces in this particular case the determinan- 
tal construction of the Jack polynomials found by Lapointe, Lascoux, and Morse 
|TXM2l . 

In dealing with the type A root system, it is more convenient to work with parti- 
tions rather than with the weights themselves. Let be the set of partitions with 
at most N parts, i.e., the set of weakly decreasing TV-tuples with components given 
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by nonnegative integers. For A,/i G An the dominance order on these partitions is 
defined as 

k 

A h [i <^ |A| = \n\ and ^(Aj - Hj) > for fc = 1, . . . , TV - 1. (5.8) 

j=i 

We write A for the orthogonal projection of A G An onto the hyperplane E C M. N 
perpendicular to the vector (1,1, . . . , 1): 

A= (A 1 ,...,A JV )-^(1,...,1). (5.9) 

The cone of dominant weights associated to the root system An-i is now given by 
the projection of An onto the hyperplane E, i.e. V\ = {A | A G A^r}, equipped with 
a partial order induced by the dominance ordering of the partitions in Eq. Ij5.8|l . 
Specifically, for given A G Ajy all dominant weights smaller or equal to A G V\ are 
given by 

Vf A = {ft | /i G A N A /i ^ A}. (5.10) 

The projection A — » A <|5.9|) has a nontrivial kernel of the form (1,1,..., 1)N. The 
set in Eq. (|5.1U|) . however, is clearly independent of the particular choice for the 
partition A projecting onto the dominant weight A. 

The Weyl group acts transitively on the root system Ajv-i, as all roots have the 
same length. Thus, the value of the root multiplicity parameter g a does not depend 
on a, viz. g a — g for all a G R. Given a partition A G Ajy, let us define for fj, ^ A 

N 

4 = E H + 5 + 1 - 2j)), (5.11a) 

and for i> ~< fi ^. \ 

!2g (mi - m 2 )N v {n 1 , n 2 ) if(J,Qv= ({mi, m 2 }, {ni, n 2 }) 
with mi — ni = ri2 — mi > 0, (5.11b) 
otherwise, 



where 

\C , " ( 2 ril) ) if |m = |n 2 |, V ' 

and r) v (n) denotes the multiplicity counter defined in Eq. I|5.7bfl . The super- 
diagonal ~ e f (A ^ ^) an d the lower triangular block d^ (y -< /t ^ A) of the 
Hessenberg matrix in Theorem 15.41 become for the .Ajv-i-type Heckman-Opdam 



A. 



polynomial 



4-4=4-4 and 4> = <- ( 5 - 13 ) 



respectively. (Notice in this connection that the expressions — and d^ v on 
the r.h.s. are invariant with respect to the additive action of (1, 1, ... , 1)N on An-) 
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5.3.2. The case Dm. The cone of dominant weights "P^ consists of the iV-tuples 
A = (Ai, . . . , Aj\r) with parts \j that are all integers or all half-integers subject to 
the ordering 

Ai > A 2 > ••• > Aiv-i > |Aiv|. (5.14) 
The partial order on V+ is defined as 



A t i" 



E-=i(Ai-Mi)eN 
Efci X (^ - N) + s(X N - Miv) G 2N for e = ±1 



for k = l,...,N-2, 



(5.15) 



The Weyl group again acts transitively on the root system Dn, so we have 
9a = J, Va £ R. The super-diagonal — (fx ^ A) and the lower triangular 
block dP v (y -< ^ A) of the Hessenberg matrix in Theorem 15.41 become for the 
Djv-type Heckman-Opdam polynomial of the form 

N 

e£=5>i(/^ + 2<?(iV-j)) 



(5.16a) 



,;=i 



and 



(^mi,m2;ni.n2 ^mi ,rri2',ni ,772 ^ ^m 1 .m 2 ;ni,n2 ^mi ,7772 ;ni ,ri2 ) (^1 1 ^2 ) 

if fi Q v — ({mi, m 2 }, {111,712}) with 71^(0) ^ and m 2 ^ 0, 

( t 'mi,ra ! ;iii,n2 + ^mi ,m^;ni ,772") ^v{ n li n 2) 

\i liQ v — ({mi, m 2 }, {^li "-2}) with n^O) = or m 2 = 0, 



<,aT ;A+ ^(A + ; 



n)+d A — A AC(A",n) 

' m,A-;A-,n v ' ' 



if ^ e v = ({m}, {n}) with n M (0) 7^ 0, 



m,A+;A+,n y 



i{fj,Qv= ({m}, {n}) with n^O) = 0, 
otherwise. 



Here m . ni „ 2 refers to the Ai-type matrix elements (cf. Eq. (|5.11b|) '). viz. 



(5.16b) 



mi .7712 ;ni ,n2 



2 g (mi — m 2 ) if mi — rti = ?i 2 — m 2 > 0, 
otherwise, 



(5.17) 



and A/"„(ni,n 2 ) is the same as above (cf. Eq. (|5.12|l % ). Furthermore, m stands for 
— m and A ± = (m ± n)/2. 

Remark. In the first line of d®, at most two terms can be nonzero if n 2 = 0, and 
at most one term otherwise. Similarly, in the second line, at most one term can be 
nonzero. 
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5.3.3. The case Bn- The cone of dominant weights Vg consists of the TV-tuples 
A = (Ai, . . . , Aj\r) with parts \j that are all integers or all half-integers subject to 
the ordering 

Ai > A 2 > ■• • > Aat_i > X N > 0. (5.18) 
The partial order on Vg is defined as 

fe 

A t A« ^ ^(Aj - nj) e N for k = 1, . . . , N. (5.19) 

i=l 

The i?Ar-type root system has two root lengths, so the action of the Weyl group 
splits up in two orbits. We will denote the root multiplicity parameters for the long 
and short roots by g and g s , respectively. The super-diagonal — di A) and 

the lower triangular block d^ v (v -< /.i H A) of the Hessenberg matrix in Theorem 
l5.4l become for the B^-type Heckman-Opdam polynomial pf of the form 

JV 

«£ = l>j(A*i + 2ff(JV -j)+g s ) (5.20a) 

3=1 

and 



with <ff v taken from Eq. (|5.16b|) and 



C = { n n u ♦ (5.20b) 



2g s mr\ v {n) if /i Q v = ({m} , {n}) 

with m — r. 
otherwise. 



d£? rt = <( with m - n > 0, (5.21) 



5.3.4. TTie cases Cat and BCn ■ The cone of dominant weights V^ c consists of the 
partitions A = (Ai, . . . , Xn) in An (cf. the Ajv_i-type above). The partial order 
on Vg C is the same as in the -B^-case (cf. Eq. (15.1911 ) 

k 

A h \i ^(Aj - nj) > for k = 1, . . . , iV. (5.22) 

j=i 

The -BCW-type root system has three root lengths, so the action of the Weyl 
group splits up into three orbits. We will denote the root multiplicity parameters 
for the long and short roots by gi and g s , respectively. The parameter for the 
remaining (i.e. middle) roots is g. The super-diagonal e^ c — ef c (/i ^ A) and the 
lower triangular block d^ (v ~< fx ^ A) of the Hessenberg matrix in Theorem 15.41 
become for the -BCjv-type Heckman-Opdam polynomial pf c of the form 

AT 

c 



and 



I ° = E N 0*j + 2 5 (N - j) + g s + 2gi) (5.23a) 



{ lD i jshort — I — ^long -c _ ,-j 

> + : 7\ , (5.23b) 
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where d® v and d s ^° rt are taken from Eqs. 15.16bjl and (|5.21(l . respectively, and 

{Agimr] v {n) if /j 9 v = ({m}, {n}) 
with to - n E 2N, (5.24) 
otherwise. 

Remark (i). The Cat case is obtained from the BCn case by setting g s = 0. In this 
situation one generally can reduce the size of the Hessenberg matrix, as the partial 
order on the weights for the CV root system, viz. 

\E J =i( A i-^)e2N, 

is less refined than the partial order in Eq. (|5.19|) corresponding to the BCn root 
system. More specifically, if the monomial on the I th row is not comparable to the 
leading monomial in the Cn ordering (|5.25|l . then we may eliminate (for g s = 0) 
the I th row together with the (I + l) th column from the Hessenberg matrix. (To 
keep the normalization monic, we should of course also delete the corresponding 
factors from the normalization constant £\.) 

Remark (ii). The BCn Heckman-Opdam polynomials deserve special attention as 
they are universal in the sense that the polynomials associated with the other 
classical root systems can be obtained as special cases (the types Bn, Cn and Dn 
by specialization of the parameters gi and g s , and the type An-i by selecting the 
top-degree homogeneous component). For a systematic study of the properties of 
the SCV-type Heckman-Opdam polynomials we refer to Refs. |BOIID2| and papers 
cited therein. 

Remark (iii). Example: for R = B3 and A = (2,1,0) the determinantal formula 
reads 

P2,l,0 = ((2 + 4g)(l + 2g + g s )(3 + Ag + g s )(4 - 

mo,o,o — 5 — 10g~3g s 

rai.0,0 6g s -A-6g-2g a 

Wl, 1,0 24g Ag 3 - 

"12,0,0 12 g s Ag 3 

"ii,i,i 8 <j 

"12,1,0 24g + 8g s 

This polynomial may also be interpreted as a special case of the £?C3-type Heckman- 
Opdam polynomial^, 1,0 with gi — 0. We observe in this connection that for g s = 0, 
the 1 st , the 3 rd and the 4 th row, together with the 2 nd , the A th and the 5 th column, 
may be eliminated from the Hessenberg matrix (cf. Remark (i) above). Indeed, 
the weights (0, 0, 0), (1, 1, 0) and (2, 0, 0) are not comparable to the highest weight 
(2,1,0) with respect to the CAr-type partial order in Eq. (|5.25() . (To keep our 
normalization monic, we must also delete the 2 nd , the 3 rd , and the 5 th factor from 
the normalization constant.) 

6. Macdonald Polynomials: the case t a = t 

In this section we apply the formalism of Section [3] to arrive at a determinantal 
construction of the Macdonald polynomials. Throughout this section it will be 
assumed that the root system R is reduced and that the dual root system i? v has 
a minuscule weight (thus including the types An, Bn, Cn, Dn, Eq and E7 while 



6g + 2 9s )(5 


+ Wg + 3g s )) 


-1 

X 




















-4p -g 3 








49 


-1 - 2g - g s 





2<? s 





-2-4 


Sgs 


4 g s 


12 g 
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excluding the types BC n , Eg, F4 and G^)- We will furthermore restrict to the case 
that t a = t, Va G R and — unless explicitly stated otherwise — we will consider the 
(q, t) parameters as indeterminates rather than real (or complex) numbers. 

6.1. The Macdonald operator. For x £ E, we define the g-translation in A via 
its action on the exponential basis: 

T x , q e X = q^e X (A £ V). (6.1) 

Definition ( M2 M3 ). Let n be a minuscule weight for i? v , i.e., the vector 7r G E 
is such that (n, a) £ {0, 1} for all a £ R + . The q-difference operator 

D * = ^ E ( II !- e ^) ) T ^>« (6 - 2) 

1 1 w€W a£R+ 

is called the Macdonald operator associated to the minuscule weight n. 

(The above definition of the Macdonald operator D n is not precisely the same 
as the one employed by Macdonald M2 M3 ; both definitions do coincide upon 
restriction to the space of invariant polynomials A w though.) In order to compute 
the action of on the monomial basis we will make use of the Weyl characters 
Xx,\£V: 

Xx = 5- 1 det H e w{x+p) , (6.3a) 

where p and S denote the half sum of the positive roots and the Weyl denominator 
respectively 

P = \ a > ( 6 - 3b ) 

aeB+ 

S = Y[ (e Q/2 -e- Q/2 ) = det{w) e w{p \ (6.3c) 

aeR+ w£W 

(Clearly the determinant det(u>) is equal to (— 1) £ ( K '), where £(w) represents the 
length of the (shortest) decomposition of w into a product of simple reflections.) It 
is well-known that for A £ V + one has that 

Xa = Kx^rrifj,, K\\ = 1, (6.4) 

with coefficients £ N. (In fact, the coefficients Kx^, which are also known as 
Kostka numbers, count the multiplicity of the weight /i in the irreducible represen- 
tation of the Lie algebra corresponding to the root system R with highest weight 
A.) An efficient way to compute the coefficients Kxu is through the application of 
Theorem 15 .41 with g a — 1, Va G R. However, for our purposes such a calculation is 
not necessary as we need the inverse of this basis transformation rather than Eq. 
itself (cf. Corollary E3 below). 
It is evident from the expansion in Eq. (|6.4|l that the Weyl characters {x\}\ev+ 
form a basis of A w . The following lemma provides a formula for the action of the 
Macdonald operator D v on the monomials m\ in terms of Weyl characters. For 
the root system R = An the formula in question is due to Macdonald |M1IIM2| . 
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Lemma 6.1 (Action of the Macdonald Operator). Let A G V + . Then one has that 
D n m x =t^ ( E t {T ' p) Q {T ' u) )xv 

Proof. Our starting point is the Weyl denominator formula in the form 
e- p Y[ (e a -l)= J2 det{w) e w(p) . 

aeR+ w£W 

By acting on both sides with the t-translator T Vi t we obtain 

t -(*, P ) e -P Y[ [&-^e a - 1) = det{w)& Mp ^e w{p) . 

a£R+ w£W 

Division of the latter identity by the former gives rise to the following expansion 
for the coefficients of the Macdonald operator 

JJ 1 - tfo°>e° = g _ H{Vtp) ^ detH ^,-(P)) e -(P). 

Substitution of this expansion in the definition of D n (taking into account the anti- 
symmetry of the Weyl denominator w(5) = det(i<;) 6), and acting on the exponential 
e v yields: 

D.e" = ^WT E det{wiW 2 )t^' W2 ^q< Wl ^' v '>e v+WlW2 ^ 

WI.W2GW 

wew, t£W(-k) 

Summation over v in the orbit W(A) then entails: 

D^m x = 8- x &' p) det{w)t {T - w(p)) q {T - u) e u+w{p) 

= 5- l &' p) det{w)t {TMp)) q {TMv)) e w{v+p) 



which completes the proof. □ 

For A G P, let w\ be the unique shortest Weyl group element such that w\(X) G 
T J+ . Then it follows from the definition of the Weyl characters that for v G V 

Xv = | det (^+p) X( w „+ p {v+p)-p) if \W V+P \ = 1, , g 5 ^ 

\0 if \W V+P \ > 1. 

(Notice in this connection that — in view of Corollarv lB.3l in Appendix[B] — the stabi- 
lizer of a weight is nontrivial if and only if there exist a root a G R + perpendicular 
to it, i.e., if and only if there exists a reflection r a , a G R + stabilizing the weight in 
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question.) Hence, to find for A, p G T J+ the multiplicity of Xp. m m\, we have to 
collect all terms in the formula of Lemma ffi. II corresponding to weights v G W(X) 
such that Wy^piv + p) — p = p, or equivalently, v — w~+ p (p + p) — p. Clearly the 
action of is triangular |M2I IMS) , since 

p = w u+p {v + p) - p = w v+p (u) -(p- w v+p (p)) < w v+p (v) < w v (v) = A 

(where in the two last steps we used the fact that any dominant weight A is maximal 
in its Weyl orbit, i.e., w{\) -< A for all w G W [Huj ). We thus arrive at the following 
explicit triangular expansion of D w m\ in terms of Xp . 

Proposition 6.2 (Triangular Expansion). Let A G V + . We have that 

pev+.p-<\ 

with 

£A = t {T ' p) q {T ' x) ; 

reW(Tr) 

h\fi = ^ det(w p+u ) e v . 

vew(\)r\(W{p.+ P )-p) 

For t = 1 the Macdonald operator trivializes to YlreWM Tr,q, which acts 
diagonally on m\ through multiplication by the eigenvalue X)rew(7r) <Z^ T '^ ■ The 
formula of Lemma 16 . II reduces in this case (and upon division by the eigenvalue) to 
the following well-known relation between the symmetric monomials and the Weyl 
characters: 

mx = Xn- ( 6 - 6 ) 

p£W(\) 

In the same way, one recovers from Proposition 16 . 21 the inverse of the expansion in 
Eq. (IOI) . 

Corollary 6.3 (Inverse Kostka Numbers). Let A G V + . The expansion of the 
symmetric monomial m\ in terms of Weyl characters is given by 

m\=X\+ J! a ^Xp, 
pev+,p^x 

with 

dA/i = ^ det(w p+u ). 

uew{\)r\(W{p.+p)- P ) 

The following proposition guarantees that the Macdonald operator is regular. 

Proposition 6.4 (Regularity). The Macdonald operator is regular in the sense 
that 

V\peV + : p^X=^e^(q,t)^e x (q,t) 
(as (analytic) functions of the indeterminates q and t). 
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Proof. After setting t — q 9 and q = exp(z), we get 
t-^P) ex = exp(z(r,A + .gp}) 

t£W(tt) 

= \W(7:)\+c 7T z 2 (X + gp 1 X + gp} + 0{z 3 ), 

with > 0. (In the last step we employed the fact that the W-invariant linear 
form X)r6W(7r)( T ' x ) vanishes and that the V^-invariant positive quadratic form 
ErewtTr) ( T > x ) mus t be proportional to because the representation of the 

Weyl group on E is irreducible and unitary). When g is positive, one has that 
(p + gp, p + gp) < (A + gp, X + gp) for all dominant weights p, X with p -< A in view 
of Proposition 15.31 It thus follows that for comparable dominant weights p ^ X 
the corresponding eigenvalues e^{q,t) and e\{q,t) cannot be equal as (analytic) 
functions of the indeterminates q and t. □ 

6.2. Determinantal construction. We will now apply the determinantal formal- 
ism of Section|31to construct the eigenbasis of D v . To this end we pick for the second 
basis {s\}xev+ the basis of Weyl characters {x\}\ev+- Specifically, by plugging 
in the eigenvalues e\ and off-diagonal matrix elements 6a m from Proposition 16.21 
together with the inverse Kostka numbers from Corollary 16.31 the formulas of 
Theorem l3.1l and the Corollaries l3 . 21 and !3 . 31 give rise to the desired eigenbasis of the 
corresponding Macdonald operator D^. For parameters such that < q, t < 1, this 
eigenbasis coincides with the Macdonald polynomials defined in Section 0] through 
the (partial) Gram-Schmidt process |M21 IM3| . We thus end up with the follow- 
ing determinantal construction of the Macdonald polynomials (not necessarily with 
0< q,t < 1). 

Theorem 6.5 (Determinantal Construction). For X G V + , let 

px = m\ + 22 C V m M 

denote the (monic) Macdonald polynomial with t a = t, Va G R. Then we have — 
upon setting for p, v G V + 

d^v = 2J det(w p+K ) (e K - e x ) 

KeW(n)n(W(u+p)- P ) 

(so d w =e ll - ex)— that: 

i) the polynomial px is represented explicitly by the determinantal formula in 
Theorem \S.l\ 

ii) the coefficients c x ^ of its monomial expansion are generated by the linear 
recurrence in Corollaru \3.''A 

hi) the expansion coefficients cxfi are given in closed form by the formula in 
Corollarv Wl\ 
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Remark (i). To determine the matrix elements d^, it is not efficient to compute the 
intersection W(p) D (W(v + p) — p) for each p, v € 7- ,+ such that v < p <\. Indeed, 
because the matrices at issue are sparse, a better strategy is to construct the matrix 
row by row. For this purpose one first determines for each dominant weight p -< A 
the set Aj, = {k G P | k £ p + W(p), \Wa\ = 1} (i.e. all regular points of the 
translated Weyl orbit p+W(p)). Weyl-permuting the weights in A M to the dominant 
cone and translating over — p produces all the nonzero contributions to the row p. 
Specifically, the nonzero matrix elements on the row corresponding to p occur in the 
columns corresponding to v from the set {y G V + \ v = Wf-(K) — p, k 6 A^}. The 
matrix elements in question are built of contributions of the form det(u>K)(es-p— ex), 
k e A M . 

Remark (ii) . It is well-known that for t — q 9 and q — * 1 the Macdonald polynomial 
Px tends to the corresponding Hcckman-Opdam polynomial (with g a — g, Va € R) 
|M3j . To perform this limit at the level of the above determinantal construction, 
it suffices to determine the asymptotics of the eigenvalues e M for q — > 1. The 
asymptotics in question is given by (cf. the proof of Proposition l6.4|) 

e^r^ = \W(ir)\ + c n (p + gp, p + gp){q - if + 0((q - l) 3 ), (6.7) 

where is a positive constant that does not depend on p and g. Since the formulas 
of the determinantal construction for the Macdonald polynomials are invariant with 
respect to an affine rcscaling of the spectrum of the form e M — > ae M + b (with a^O), 
we only pick up the second-order term of the asymptotics in Eq. I|6.7|l when sending 
q to 1. The upshot is that by replacing by {p+gp, p+gp) in Theorem l6.5l we wind 
up with an alternative determinantal formula for the Hcckman-Opdam polynomials 
(with g a — <?, Va £ R, and with R such that i? v has a minuscule weight). From 
a practical point of view the formulas coming from Theorem 16.51 are less efficient 
than those of Theorem 15.41 however, as the action of the Macdonald operator 
expanded in Weyl characters tends to be much less sparse than the action of the 
hypergeometric differential operator expanded in monomials. As a consequence, the 
matrices entering the determinantal formulas for the Heckman-Opdam polynomials 
coming from Theorem 16.51 are much less sparse than those of the type given by 
Theorem ICTI 

Remark (iii). The recurrence relation in Part ii) of Theorem 16 . 51 reads concretely 

( £ A - e xi i-i)) C AA( {-i) = 
n 

$3 det(w p+K ) (e K - ex) c xx(k} . 

k=e K ,ew(x( k ))n(W(x( e - 1 )+p)- P ) 

This relation should be regarded as a symmetrized Freudcnthal type recurrence for 
the coefficients in the monomial expansion of the Macdonald polynomials. For e M = 
(p+gp, p+gp), this recurrence degenerates to a recurrence for the coefficients of the 
Heckman-Opdam polynomials with g a — g, Va G R (cf. Remark (ii) above). The 
recurrence in question is different from the previous recurrence for the Heckman- 
Opdam polynomials originating from the hypergeometric differential operator (cf. 
Remark (v) at the end of Section IH~5|> . In particular, for g — 1 this gives rise to 
an alternative system of symmetrized Freudenthal type recurrence relations for the 
weight multiplicities of characters of simple Lie groups. 
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6.3. Tables for the classical root systems. We now provide tables for the con- 
struction of the Macdonald polynomials associated with the classical root systems. 
The minimum information needed for invoking Theorem 16 . 51 consists of i. the cone 
of dominant weights and its partial order, ii. the half-sum of the positive roots 
p, Hi. the action of the Weyl group, and iv. the eigenvalues e p . Below we list 
items ii.-iv for R = Ajv-i, -Bat, Cn and Dn- For item i the reader is referred to 
Subsection 15. 31 

6.3.1. The case An~i- The An~i -type Macdonald polynomials amount (in essence) 
to the Macdonald symmetric functions of Ref. |M1| . Theorem 16.51 reproduces in 
this case the determinantal construction of the Macdonald symmetric functions due 
to Lapointe, Lascoux, and Morse |LLMlj . 

We will again formulate the construction in terms of partitions with at most N 
parts, by adding a trivial center to the weight lattice (cf . Subsection I5.3fl . The 
Weyl group W is given by the permutation group of N letters Yin- A Weyl group 
element w = a G £_/v acts on a partition A G An by rearranging its parts 

(T(Ai, . . . , Ajv) = (Ao-j , . . . , Xa N )- (6.8) 

To construct the Macdonald symmetric function associated to a partition A G 
An, one employs Theorem 16 . 51 with V + = An endowed with the dominance order 
in Eq. Ij5.8jl . and 

N 

4 = E^V"' (6-9a) 

3=1 

p A = (JV-1,JV-2,...,1,0). (6.9b) 

For a rearrangement n of a partition p G An, the sign det(w p + K ) is given by the 
signature of the shortest permutation a p + K rearranging p + k such that its parts 
become weakly decreasing (i.e. det(w p+K ) = (— 1) £ ( ct p+«), where £(a p+K ) denotes 
the number of transpositions of the permutation). 

Projection of the resulting Macdonald symmetric function onto the space of ho- 
mogeneous functions of degree zero (i.e., replacing m p by m p in the monomial ex- 
pansion) entails the ^Ar^-type Macdonald polynomial associated to the weight 

A (ET31) . 

6.3.2. The case Bn- The Weyl group is the semi-direct product of Eat and the N- 
fold product of Z 2 = Z/(2Z), i.e., W = Ejv K 1^ . A Weyl group element w = (a, e) 
acts on a weight A G V\ as 

iu(Ai, . . . ,Ajv) = (£iA CTl , . . . ,£jvA CTjv .), (6.10) 

with Sj G {1, -1} for j = 1, . . . , N. 

To construct the B/v-type Macdonald polynomial associated to a weight 
A G V+, one employs Theorem 16 . 51 with 

4 = ^{t^q^+t^q-^), (6.11a) 

3=1 

pb = (JV-i,iV-|,...,|,i). (6.11b) 
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For a weight K £ W{p) with p € Vg, the sign det(w p + K ) is given by the signature 
of the shortest permutation rearranging p + k such that the absolute values of its 
parts become weakly decreasing, multiplied by (— l)" e where 

n e = \{j = l,...,N\p j + K j <Q}\. (6.12) 

6.3.3. The case Cn- The Weyl group and its action on a weight A € Vq are the 
same as in the Bjv-case. 

To construct the CV-type Macdonald polynomial p^ associated to a weight A E 
7- + , one employs Theorem 16.51 with 



N 

C 



n(i W+1 -V> /2 + <T^ /2 ), (6.13a) 



pc = (N,N -1,..., 2,1). (6.13b) 

The sign det(u> p + K ), for k £ W(p) with p £ V+, is computed in the same way as 
in the B^-case. 

6.3.4. The case Dn- The Weyl group is given by W — £jv K ", and the action 
of w = (a, e) £ W on a weight A £ is given by Eq. (|6.1U|) with Sj £ {1, -1} for 
j = 1, . . . ,N such that £i£2 4 • ■ £tv = 1- 

To construct the D^-type Macdonald polynomial p® associated to a weight 
A G , one employs Theorem 16 . 51 with 



N 

^ = ' V • : '/ >• (6.14a) 



/9 = (N-1,N-2,...,1,0). (6.14b) 

For a weight k £ with /x £ Vp, the sign det(w p+K ) is given by the signature 

of the shortest permutation rearranging p + k such that the absolute values of its 
parts become weakly decreasing. 

Remark (i). Example: for R = D% and A = (2,1,0) the determinantal formula 
reads 



ft D D w D D w D D 

£>2,1,0 = ([ e 2,l,0 - ei,l,xK e 2,l,0 ~ e X,l,-lA e 2,l,0 ~ e i,o,oJJ 



X 



mi >0 ,o £2,1,0 - £°o,e 

"11,1,-1 ^^2,1,0 + e i,-i,i e ?i,o — e M,-i 

m M,i ^ £2,1,0 + e^-i,-i £ 2^i,o — e ?,l,i 

—, -D , -D o-Ci-D i-C D,D|D 

"12,1,0 — e l,-2,0 + e -l,2,0 ~~ * e 2, 1,0 T El, 0,-2 T to, 2,-1 ~~ z e 2,l,0 + e l,0,2 T e 0,2,l 

with €° i m2 m3 = (t 4 q mi + q~ mi ) + (t 3 q m2 +tq- m *) + (t 2 q ms + t 2 q~ m3 ). 

Remark (ii). For the root systems Bn and Cn the minuscule weights 7r of the dual 
root systems are unique. Specifically, for R — Bn the minuscule weight is given by 
lowest fundamental weight w\ of i? v (= Cn), and for R — Cn it is given by the 
highest fundamental weight ujn of R v (— Bn)- In the cases of the root systems 
An -i and Dn the above formulas for the eigenvalues e M correspond to picking for 
the minuscule weight 7r the lowest fundamental weight u)% of R v (= R). However, 



24 



J.F. VAN DIEJEN, L. LAPOINTE, AND J. MORSE 



in these cases there exist actually several alternative possibilities for the choice of 
the minuscule weight 7r. For R = An-i wc could work with each of the eigenvalues 

<r= E IK^, r = l,...,N-l, 

JC{l,...,N}jGJ 
\J\=r 

corresponding to the fundamental weights w\, . . . , wjv-i, respectively. For R = 
we could alternatively work with the eigenvalues 



&v-i = (*£+-e£-)/2 or e° N = (e° + + e°_)/2, 



where 



N 



AT 



i=i 

corresponding to the fundamental (spin) weights lon-i and ujn, respectively. From 
a computational point of view, however, in this last case it is more efficient to 
work simply with the eigenvalues + , corresponding to the linear combination of 
Macdonald operators D = D UN1 + D U{i . For instance, in the example of Remark 
(i) above, this amounts to replacing the eigenvalues m2>m3 by e^ mj m3i+ = 

g -(mi+m2+m 3 )/2^ g mi + l}(tq m2 + l)(q nl3 + 1). 

7. Macdonald Polynomials: the case of general t a 

In this section we will briefly indicate how to generalize the results of the previous 
section to the case of Macdonald polynomials with general parameters t a such that 
tw(a) — ta for all w £ W. We will keep the restriction that our root system R is 
reduced and that the dual root system i? v has a minuscule weight n. 

For general W-invariant i Q -parameters the Macdonald operator becomes |M3| 

1 / 1 — t^' a "> e w i a ) X 



1 1 - 



It is convenient to reparameterize the i Q as 

t a = q 3a 

(with ^(q.) = g a , Vw € W). The action of A on the monomial basis can be 
written as M3, Section 5] 

A,** =«<*•>•> E (-l) |x| £ ^^ (XC) ^ (X)) X^ 2pm , (7-2) 
xcr+ uew(x) 

where X c = R+\X and 

(so p = p(R + ) and p 9 = p g (R + )). Bringing the action in Eq. (|7.2|l to triangular 
form gives 

At ™A = £A YA + E 

(7.3) 
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with 

ex = q<"'") ]T q (T ' X+p3 K 

t£W(tt) 

£ det(«, w+p(Jc . ) - pW )ff (,r ^ <jr) -^ OT >. 
i»EW(\)n(w( f i+p)+p(x)-p(x c )) 
XcR + 

We thus wind up with the following determinantal construction of the Macdonald 
polynomials for general VF-invariant parameters. 

Theorem 7.1 (Determinantal Construction). For A G V + , let 

P\=m\+ ^ c ^ "V 

denote the (monic) Macdonald polynomial with t a — q 9a . Then we have — upon 
setting for /j,, v £ V + 

t£W(tt) 

d^iv — bpv e\a^ij^ 

with 

eipu = ^ det(w p+K ), 

= 

£ dat(«; M . p(x .,_ p(x) )ff<-.-+''.( Jc0 )-^W> 
Kew(n)n(W(v+p)+p(X)- P (x e )) 

XcB + 

(so = ep, - e\)— that: 

i) the polynomial px is represented explicitly by the determinantal formula in 
Theorem \H.l\ 

ii) the coefficients cxp, of its monomial expansion are generated by the linear 
recurrence in Corollaru \3.''A 

hi) the expansion coefficients cxp are given in closed form by the formula in 
Corollarv Wlk 

Remark (i) . From a computational standpoint the formulas of Theorem 17.11 are 
much less effective than the determinantal constructions for the (q, t) Macdonald 
polynomials fTheorem l6.5|l and (especially) for the Heckman-Opdam polynomials 
(Theorem 15. 411 . This is because the action of the general Macdonald operator on 
the monomial basis (cf. Eq. I|7.2|l ) is much less sparse and the matrix elements 
are moreover much more complex than in these two previous cases. This renders 
Theorem 17. II presumably only of limited practical value. 

Remark (ii). The most general class of Macdonald polynomials admits a richer 
parameter structure connected with admissible pairs of root systems (R, S) |M3j . 
(From this perspective the polynomials studied here are of the type (R, R).) Since 
Macdonald in fact gives the action of the Macdonald operator on the monomial 
basis for general admissible pairs, it is not difficult to generalize Theorem 17.11 also 
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to this context (at the expense of having to introduce a more elaborate notational 
apparatus). 



Appendix A. Determinants of Hessenberg Matrices 

In this appendix we recall a classic recursive method for the efficient evaluation 
of the determinant of a Hessenberg matrix . This recursive method was used in 
Section|2|for the explicit evaluation of our determinantal formula for the eigenbasis 
of the regular triangular operators in A. . 

Lemma ( W ). Let \D\ be the Hessenberg determinant 

mi —di2 • • ■ 



\D\ 



m 2 



l 2,3 



dn-1,2 
d n .2 



dn-1,3 
dn,3 






dn ,n 



with nonzero elements on the super- diagonal: dj—ij (1 < j < n). Then the 
expansion of \D\ with respect to the first column is of the form 



\D\ = ^ c i m i> 



with cVi 



di, 2C?2,3 ■ • • d n -i„ and 



Ci-l 



3=i 



Cj dj t t 



(K£<n). 



Proof. ( LLM2|) We denote the columns of our Hessenberg matrix by the n-dimen- 
sional vectors m and d^ fc ^ (1 < k < n), respectively. Let c — (ci, Ca, . . . , c„) 
be a nonzero vector perpendicular to the hyperplane spanned by the (n — 1) lin- 
ear independent columns d^ 2 ), . . . ,d^ n '. Because the value of the determinant is 
equal to the (hyper)volume of the polygon determined by df 2 v..,dW, multi- 
plied by the height of m in the perpendicular direction c (possibly up to a sign), 
we conclude that \D\ must be proportional to the scalar product of c and the 
first column m, i.e. \D\ ~ (c, m) = ^™ =1 qto£. Furthermore, since the cofac- 
tor of m n in \D\ equals the product of the elements on the super-diagonal, one 
sees that in fact |D| = Y^e=i c £ m e provided the normalization of c is fixed such 
that c„ = di 2^2,3 ••• ^n-i,n- The lemma now follows from the observation that 
the requirement that c be orthogonal to the columns d^ 2 \ . . . , d^") translates itself 
directly into the stated recurrence relations for the components c,, j = 1, . . . , n. □ 



Appendix B. Counting the Orders of Stabilizer Subgroups 

To build the determinantal formula for the Heckman-Opdam polynomials in 
Theorem 15 .41 one frequently needs to compute the orders of stabilizer subgroups of 
the Weyl group. In this appendix we include a short proof of a useful formula for 
the orders of these stabilizers that can be found e.g. in Ref. [MSI Section 12]. 
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Proposition B.l (Orbit Size). Let A e V . Then the size of the Weyl orbit through 
A is given by 

„ (p.a v ) + 1 + U% 

\W(X)\= n 



<A,a v }#0 

where p—\ J2aeR+ a ' an d 

1 otherwise. 

Proof. Let us first assume that A is dominant. The orbit size can be obtained in this 
case from the following evaluation formula for the Heckman-Opdam polynomial p\ 
at x = (i.e. at the identity element of the torus T = £/(2ttQ v )) 

rrn TT I^9' aV ) +9" + 55f ]<A,aV) 

where [a] m = a(a + 1) • • • (a + m — 1) (with the convention that [a]o = 1), and by 
definition g Q = if a $ R. Indeed, setting g a = g for all aeii, and performing 
the limit g — > 0, readily entails the formula of the proposition. (Here one uses 
that lim gcv _^oPA = m A an< i that m^(0) = |W(A)|.) The extension to non-dominant 
weights A is immediate (cf. also the two corollaries below). □ 

By picking A regular (for instance strongly dominant), one gets a formula for the 
order of the Weyl group. 

Corollary B.2 (Order of the Weyl Group). The order of the Weyl group is given 
by 

(p,a v ) + l + U % 



\w\= n 



aeR + V" / 1 2 

Dividing the order of the Weyl group by the size of the orbit produces the order 
of the stabilizer subgroup. 

Corollary B.3 (Order of the Stabilizer). Let A e P. The order of the stabilizer of 
A is given by 

\w x \= TT (p ' aV) + 1 + |(5f 

<A,Q v )=0 

Remark (i) . It is clear that the formulas of Proposition IB. II and Corollary IB. 31 in 
fact serve to compute the sizes of the orbit W(x) and the stabilizer W x for any 
vector x 6 E (i.e. not just weight vectors). 

Remark (ii). For a reduced root system R the above formulas simplify somewhat 
as 8<* = in this case. 

2 
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